| 学校名 | 生徒数 | 平均点 | 抽出確率 | 標本ウェイト |
|---|---|---|---|---|
| A校 | 150人 | 50点 | \(\frac{20}{150}\) | \(\frac{150}{20}=7.5\) |
| B校 | 100人 | 60点 | \(\frac{20}{100}\) | \(\frac{100}{20}=5.0\) |
| C校 | 50人 | 70点 | \(\frac{20}{50}\) | \(\frac{50}{20}=2.5\) |
| D校 | 20人 | 80点 | \(\frac{20}{20}\) | \(\frac{20}{20}=1.0\) |
| 全体 | 320人 | ?点 |
★全体の平均点をどうやって計算する?
⇒ 母集団の推定が(超)面倒くさい
# データの読み込み
url <- "https://raw.githubusercontent.com/kawa5902/LSAdata/refs/heads/main/pisa2012stuJPN.csv"
jpn2012 <- read.csv(url)
# ESCSの欠測をNAに変換
jpn2012$ESCS[jpn2012$ESCS == 9999] <- NA
# ESCSの平均値
intsvy::pisa.mean(variable = "ESCS", data = jpn2012) Freq Mean s.e. SD s.e
1 6185 -0.07151303 0.01524119 0.7036702 0.007597166
ST03Q01 Freq Percentage Std.err.
1 1 521 8.20 0.36
2 2 456 7.27 0.38
3 3 501 7.89 0.38
4 4 532 8.30 0.36
5 5 533 8.43 0.34
6 6 521 8.07 0.30
7 7 574 9.00 0.40
8 8 560 8.94 0.36
9 9 541 8.56 0.43
10 10 547 8.64 0.32
11 11 510 7.98 0.38
12 12 555 8.72 0.38
ST04Q01 Freq Mean s.e. SD s.e
1 1 2974 -0.07228828 0.02079859 0.7075644 0.010526624
2 2 3211 -0.07079883 0.01925711 0.7001372 0.008863656
ST04Q01 ST03Q01 Freq Percentage Std.err.
1 1 1 263 8.73 0.51
2 1 2 225 7.60 0.51
3 1 3 226 7.48 0.49
4 1 4 249 8.11 0.48
5 1 5 268 9.00 0.49
6 1 6 247 7.97 0.42
7 1 7 263 8.57 0.58
8 1 8 264 8.78 0.49
9 1 9 249 8.40 0.62
10 1 10 275 9.16 0.50
11 1 11 226 7.43 0.47
12 1 12 266 8.76 0.51
13 2 1 258 7.73 0.49
14 2 2 231 6.97 0.46
15 2 3 275 8.25 0.49
16 2 4 283 8.47 0.46
17 2 5 265 7.92 0.48
18 2 6 274 8.15 0.45
19 2 7 311 9.38 0.51
20 2 8 296 9.08 0.50
21 2 9 292 8.70 0.50
22 2 10 272 8.17 0.51
23 2 11 284 8.49 0.49
24 2 12 289 8.69 0.52
PV1READ Rho PV1READ s.e. PV1MATH Rho PV1MATH s.e. ESCS Rho ESCS s.e.
PV1READ 1.000 0.000 0.861 0.006 0.289 0.024
PV1MATH 0.861 0.006 1.000 0.000 0.320 0.025
ESCS 0.289 0.024 0.320 0.025 1.000 0.000
attr(,"class")
[1] "intsvy.rho" "matrix" "array"
Freq Mean s.e. SD s.e
1 6351 536.41 3.59 93.52 2.19
# PVsを使った平均値の推定 by 男女別
intsvy::pisa.mean.pv(
pvlabel = paste0("PV", 1:5, "MATH"), by = "ST04Q01", data = jpn2012
) ST04Q01 Freq Mean s.e. SD s.e
1 1 3021 527.01 3.59 88.06 2.63
2 2 3330 544.88 4.62 97.41 2.65
Estimate Std. Error t value
(Intercept) 543.78 3.32 163.92
ESCS 39.24 3.69 10.63
R-squared 0.08 0.01 5.95
# PVsを使った回帰分析
intsvy::pisa.reg.pv(
x = "ESCS", pvlabel = paste0("PV", 1:5, "READ"), data = jpn2012
) Estimate Std. Error t value
(Intercept) 543.59 3.34 162.58
ESCS 37.99 3.91 9.71
R-squared 0.08 0.01 5.42
# PVsを使った回帰分析(重回帰)
intsvy::pisa.reg.pv(
x = c("ST04Q01", "ESCS"), pvlabel = paste0("PV", 1:5, "READ"), data = jpn2012
) Estimate Std. Error t value
(Intercept) 577.85 5.10 113.30
ST04Q01 -22.53 3.54 -6.36
ESCS 38.00 3.94 9.65
R-squared 0.09 0.01 6.51
1